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Abstract 

We compute the electromagnetic response and corresponding forces between two silver nanowires. 
The wires are illuminated by a plane wave which has the electric field vector perpendicular to the 
axis of the wires, insuring that plasmonic resonances can be excited. We consider a nontrivial square 
cross section geometry that has dimensions on the order of 0.1 A, where A is the wavelength of the 
incident electromagnetic field. We find that due to the plasmonic resonance, there occurs great 
enhancement of the direct and mutual electromagnetic forces that are exerted on the nanowires. 
The Lippman-Schwinger volume integral equation is implemented to obtain solutions to Maxwell's 
equations for various A and separation distances between wires. The forces are computed using 
Maxwell's stress tensor and numerical results are shown for both on and off resonant conditions. 
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I. INTRODUCTION 



The phenomenon of surface plasma resonance excitation (surface plasmon) was first 
predicted^ as an explanation of observed energy loss spectra of electron beams penetrating 
thin metal films. Also, at optical wavelengths the reflectance of metal (Al, Ag, Au) mirrors 
with surface roughness is significantly affected by surface plasmon excitation 2 . When con- 
sidering scattering of plane waves at wavelength A by particles of dimension << A, Rayleigh 
scattering theory can often be applied. In this case the shape of the particle is not important 
and the scattering cross section is typically quite small. However more recently, the study 
of electromagnetic enhancement and field localization by nano-sized metallic structures via 
plasmonic resonance has been the focus of research activity. In this case, when plasmonic 
excitation is involved, the shape and material details of the structure become important 
with respect to conditions for resonance and this greatly increases the scattering of elec- 
tromagnetic waves. When two or more nanosized structures are illuminated at a resonant 
wavelength, the electromagnetic coupling and associated charge redistribution that occurs 
can result in a significant mutual coupling force. Plasmon resonances lead to extremely large 
localized fields at specific wavelengths in the the vicinity of nanoparticles, thereby resulting 
in a large scattering cross section^ and mutual forces between nanoparticles. The localized 
fields also play an important role in the surface enhanced Raman scattering. This enhance- 
ment can be large enough to enable the detection of a single molecule. A primary driving 
force for the recent focus on plasmon resonance of nanoparticles is the application of opti- 
cal properties of these structures in novel optical devices^ and biosensors^.. It has been 
noted that nanowires with different cross sectional shapes such as, cylindrical, triangular 
and elliptical exhibit multiple resonances depending on their particular shape. For example, 
the field enhancement in the vicinity of a 20 nm triangular particle can exceed 400 times 
the incident field amplitude, while it is only 10 for a cylindrical particle of the same size£. 
In the past few years, several researchers have investigated the resonance effects of particle 
sizes larger than 100 nm. Very recently, attention has been focused on resonance effects of 
nanoparticles in the 20 to 50 nm ranged. 

A primary difficulty in dealing with particle size smaller than 100 nm, especially during 
resonance conditions, is the increasing computational complexity due to rapid variations 
of the field over short distances. For structures that deviate from the simpler cylinder 
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and sphere geometries, analytical approaches are not available, and one must rely on an 
effective numerical approach that correctly incorporates the necessary boundary conditions. 
An efficient numerical procedure for the computation of fields at the plasmon resonance 
wavelength of nanoparticles of arbitrary cross section is afforded by the Lippman-Schwinger 
integral equation^. This numerical formalism is a volume integral equation utilizing the 
Green's tensor for 2D or 3D geometries. In this paper, we calculate the electromagnetic 
forces exerted on two nanowires with square cross section. The side dimensions of the 
nanowires are 30 nm and the wavelength of the incident field varies from 300-516 nm. 

II. METHOD 

We first outline the volume integral equations, which are appropriate for the investiga- 
tion of nanosystems, and in this case electromagnetic resonances in 2D coupled nanowire 
structures. The system is assumed infinite in the z direction, and all spatial variation occurs 
in the x-y plane. The physically relevant quantities that we shall focus on are the electric 
(E) and magnetic (H) fields and the nature of net forces exerted on a given wire under 
conditions of plasmonic resonance. In the absence of magnetic media (// = 1), the basic 
integral equations are written as^, 

E(r) = E inc (r) + k 2 [ dVG e (r, r') ■ 5e(r')E(r'), (la) 
Js 

H(r) = H inc (r) - ik [ dVG m (r, r') ■ 5e(r')E(r'), (lb) 

Js 

where k = u>/c, r = (x, y), and the integration is over the cross-sections 5* of the nanowires, 
which are embedded in vacuum. The free-space electric dyadic Green's function, G e (r,r'), 
is a second-rank tensor and its explicit form has been given elsewhere^. The magnetic 
dyadic Green's function, G m (r,r'), is calculated using G m (r,r') = V x [I go(r, r')] , where 
go(r, r') = (i/4)/f (^oP) exp(ik z z). Since the incident field propagates solely in the x—y plane 
{k z = 0), this yields the following nonzero components: G™ = — z/(4p)fc (y — y')Hi(kop) 
and = i/(4p)k (x — x^H^kop), where p = |r — r'|. The other components are found 
straightforwardly using the antisymmetric relation G™j = —G™i- The scalar permittivity 
contrast is mapped according to <5e(r) = e — 1, when r is entirely within the nanowire with 
permittivity e, and otherwise. We then cast the integrand of ()la|) into the form of a 
linear matrix equation system by discretizing the wire surface as a fine grid for numerical 
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integration. Once the electric field is found within the scattering medium, E(r) is then 
inserted into (jlbj) to calculate H(r). As a check, we calculated the H field directly from 
the E field using the Maxwell-Faraday law and found consistency with the results obtained 
using Eq. (JTEJ). 

To calculate the net force acting on a particular wire, we begin with Maxwell's stress 
energy tensor, T, which has components,— 

T af3 = i- [E a Dp + H a B p - I(E • D + B • H)^] . (2) 

The net force on an object, enclosed by an area A is then given by F = {§ A {T ■ n) dA) 
where n is the outward normal to the surface A enclosing a given nanowire, and () denotes 
the time average. For the geometry under consideration, and with E = (E x ,E y ), only the 
T xx ,T xy , T yx and T yy components of the stress tensor contribute to F. 



III. NUMERICAL RESULTS 

The Lippmann-Schwinger equation (Eq. (|la|) can be cast into the standard form of a 
linear equation system Ax = b, and is solved here numerically by an efficient and stabilized 
version of the conventional BiConjugate gradient method.— The matrix A does not have 
to be stored in computer memory, allowing discretization of the scatterer over a very fine 
scale without undue limitations on the storage requirements. We model the permittivity of 
the metal nanowires using the experimental silver data in Ref. The square nanowires 
are taken to be 30 nm per side dimension and the grid count per nanowire is 80 x 80. This 
yields a grid spatial resolution of 0.375 nm x 0.375 nm. Several numerical tests confirmed 
that this resolution is quite satisfactory. 

In general for silver nanowires with small cross sections, the inverse collision time for 
electrons can increase from the additional scattering events that take place at the surface, 
and thus scattering from the boundary becomes relevant. This problem has been treated 
using alternate techniques. The effect of electron screening is expected to also reduce the 
inverse collision time, however a careful modeling of the particle surface showed^ an effective 
scattering rate that is comparable to the classical result. 

Significant attraction between two metal nanowires can be generated by Coulombic forces. 
If the gap distance is such that the plasmonic evanescent near-fields overlap, then the poten- 
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FIG. 1: (Color online) The wavelength dependence of the x-component of the time-averaged total 
force on the left wire (wire 1) for five different separations d, as indicated in the legend. The 
incident plane wave propagates downwards in the — y direction, while the field is polarized along 
the x direction, which is appropriate to excite surface plasmons. The force is always positive, 
indicating attraction between the nanowires since the adjacent wire (wire 2) feels an identical force 
that is opposite in sign (not shown). The inset illustrates the corresponding x-component of the 
dipole moment, with the solid gray curve representing the case of a single wire. 

tial for coup led- wire resonant modes exist. These coup led- wire modes can be very complex 
depending on gap distance, permittivity, wavelength, and geometric parameters. These 
modes can have a widely varying spatial distribution of the electron plasma. It then fol- 
lows that if the gap fields are enhanced because of a resonance condition, there should be 
enhanced Coulombic forces. 

To address this, we begin by examining the forces exerted on the nanowires due directly 
to the incident wave and the mutual coupling between the wires. The results naturally 
depend on the propagation direction of the incident field, given by its wavevector k = 
— ^(cos^x + sin#y). In Fig. we illustrate the x-component of the time-averaged force, 
F x , on the left wire (wire 1) as a function of the incident field wavelength, A, and for various 
gap distances d. The incident electric field vector is polarized along the x direction, and 
the wavevector is directed along — y. This particular incident field propagation direction 
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yields oppositely directed net forces on each wire that are equal in magnitude. It is evident 
that F x is positive for wire 1 over the whole A range shown, indicating attraction among the 
nanowires. It is seen that for much of the A range, smaller gaps yield greater maximum values 
of F x . All curves show a decrease and blueshift in the F x peak as d increases, except for the 
d = 20 nm curve, which shows very little apparent interaction between the wires. This is 
consistent with a coupled-wire resonance condition that is supported by the p x inset showing 
the rc-component of the dipole moment of wire 1. This quantity is obtained by integrating 
the appropriate component of the real part of the polarization, P(r) = l/47r5e(r)E(r), over 
the cross sectional area of the wire. We show here only p x since this component of the dipole 
moment encapsulates the relevant charge distribution involved in interwire coupling. In the 
inset the solid gray curve shows for reference, p x for a single isolated wire. Note that the 
d = 20 nm curve for p x is very similar to the solid gray curve, further indicating very little 
mutual interaction between the wires. For 330 nm < A < 400 nm, p x has negative slope 
and becomes more negative quite rapidly with wavelength. This is approximately the region 
over which the F x curves reach their primary peak. It is also interesting that for larger 
wavelengths, there is another set of F x peaks that are visible in the d = 3, 5, and 7 nm 
curves. These peaks become washed out at larger d, and approximately correlate with p x 
(inset) where the slopes are positive. 

The sign change in the dipole moment exhibited in Fig. ^ is easily visualized by means 
of vector plots. To this end we show in Fig. |2] the electric field vector distribution within 
the nanowires. The wires are separated by 5 nm and the top and bottom panels correspond 
to incident field wavelengths of A = 340 nm and A = 385 nm respectively. In the top 
A = 340 nm panel, the corner dipole fields are most intense on the outer left and right sides. 
There is a subsequent nonuniform, strongly position dependent field orientation within the 
nanowires. The E field on the sides nearest the gap region is relatively small compared to 
the corner fields. In the bottom A = 385 nm panel, the situation is markedly different. Here 
the gap field is dominant and different mode patterns arise. These results are consistent 
with Fig. ^ where F x at A = 340 nm is small compared to A = 385 nm. 

We now vary the incident electric field so that it is now polarized in the y direction, 
and propagates along — x. We show in Fig. |3(a)[ F x on wire 1 as a function of A. Figure 
|3(b)| depicts the same quantities for wire 2. In general, the net force on each wire is no 
longer equal in magnitude due to the uneven radiation effects directly from the incident 
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FIG. 2: Vector plots of the E field within the nanowires illustrating the local field directions. Here 
an incident plane wave is traveling downwards as in Fig. ^ The wires are separated by 5 nm. The 
top panel corresponds to A = 340 nm while the bottom panel is for A = 385 nm. The top panel 
clearly shows that the E field is most intense at the outer corners opposite the gap. Integration 
of the polarization over the left wire yields p x > as shown in Fig. ^ Conversely, in the bottom 
panel the field in proximity to the gap is most intense and integration over wire 1 yields p x < 0, 
again consistent with Fig. ^ This demonstrates the reason for the relatively large and attractive 
forces between the wires at A = 385 nm compared to A = 340 nm. 



7 



Wire 1 




388 Safl 340 368 1K0 400 420 44() 460 4H0 ?!X> 



X (inn) 

Wire 2 




3iX) .120 St) 3Q0 380 400 420 440 460 4H0 SOB 



A (nm) 

FIG. 3: (Color online) The wavelength dependence of the x-component of the time-averaged total 
force on wire 1 (a) and wire 2 (b) for five different separations d, as indicated. The coordinate 
system is identical to that shown in Fig. ^ The incident plane wave propagates with k in the — x 
direction and the electric field is polarized in the y direction. The dashed curves correspond to F x 
on a single nanowire. There are also insets that show the net dipole moment p y for each nanowire. 
For this incident field orientation p x is zero. For reference, the solid gray curves in the two insets 
shows p y for an individual nanowire. Note that in the neighborhood of A ~ 336 nm, p y < in (a) 
and p y > in (b) which yields a repulsive force between the nanowires. When A ^ 400 nm, the 
wires are largely attractive, except for d = 20 nm, where there is large enough separation that the 
interaction between wires is minimal. 
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wave. Figure p(a)| shows that for sufficiently small wire separations (d < 10 nm), the force 
on wire 1 changes sign in a continuous manner as A is varied. The force on wire 2 is mainly 
negative except for a narrow band of A in the vicinity of A ~ 336 nm (see Figure |3(b)| ), 
which correlates to when the force on wire 1 is most negative. Thus the wires experience 
a repulsion at these wavelengths and separations. For larger separations F x diminishes and 
tapers towards zero at higher wavelengths, indicating that electromagnetic coupling between 
the two wires is negligible. The plasmon resonances within these structures reveal themselves 
through prominent peaks in the force curves in both figures. The peak positions are strongly 
dependent on d, as they broaden and reduce in magnitude as d increases. Comparing the 
scales in Figs. Q and E3 we find that the overall resonance effects observed for the forces at 
small separations is reduced somewhat for this incident field direction. The insets in Fig. |3] 
show the y-component of the electric dipole moment (p y ) versus A for wires 1 and 2. For 
wire 1 and d = 5 nm, Fig. |3 (a)] shows that p y undergoes a minimum at A ~ 328 nm and then 
sharply crosses zero at A ~ 332 nm. This trend is similar for the other separations, with only 
a slight redshifting for larger d. For both wire 1 and wire 2 a clear correlation exists between 
p y and the associated force; for small enough d, the pronounced peaks or depressions in p y 
correspond to when \F X \ is largest. These features therefore become redshifted and spread- 
out with decreasing separation. As one would expect for d = 20 nm, the general features of 
p y approach that of the limiting case of a single nanowire. 

To contrast the near-field electric field patterns for differing incident wave directions, we 
display in Fig. H]the electric field amplitudes in the region of the two nanowires. The wires 
are set at a fixed distance of 5 nm and are illuminated at two different wavelengths. The 
top two panels correspond to an incident field directed downwards (as shown by the arrows) 
with A = 385 nm and A = 470 nm. The bottom two panels are for an incident electric field 
traveling towards the left. Beginning with the top left panel, at A = 385 nm there is a large 
field enhancement from the corners and within the gap between the wires. This is consistent 
with Fig. ^ where the force has a maximum at that wavelength. As the wavelength is 
increased to A = 470 nm, the system is no longer at resonance and the interaction between 
the wires has decreased. This follows from Fig. ^ where the force is reduced by nearly a 
factor of three from its peak value. A different situation is observed in the bottom panels, 
where at A = 385 nm the electric field is small between the wires, although the outer corners 
are still illuminated. At A = 470 nm, the field amplitude is intensified predominantly within 
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FIG. 4: (Color online) Electric field amplitude contours (normalized by the incident field) for two 
nanowires illuminated by a plane wave with wavevector direction indicated by the arrows. The 
separation distance is d = 5 nm. The electromagnetic response of the nanowire pair is quite different 
depending on the angle of incidence and wavelength. In the two upper plots, the electromagnetic 
response and subsequent attractive force is greatest for the shortest wavelength. Conversely, in the 
lower plots the response is greatest for the longer wavelength. These results are consistent with 
Figs. [T] and El 

the gap as the nanowires are electrically coupled. There is a clear relationship between 
the force profiles in Fig. |H] and the field amplitude plots. It is apparent that the incident 
field "sees" a different geometry based on its orientation relative to the coupled structure; 
downward propagating waves interact with an effectively larger nanowire in addition to the 
interwire gap that may support additional modes due to the wire-wire interaction. 

Next, in Fig. the force on each of the two nanowires is shown as a function of the 
separation distance d, and for an incident field directed along — x. Also shown is the cor- 
responding induced net dipole moment, p y , for each wire. There are six A considered and 
the variation of F x with d is strongly wavelength dependent. Consider first the longest 
wavelength, A = 470 nm where the separation range d = 4 — 5 nm shows a large positive 
F x in Fig. 03(a) and a correspondingly large negative F x in Fig. Efb). This indicates a res- 
onant condition resulting in a relatively strong attractive force. This is in contrast to the 
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FIG. 5: (Color online) The ^-component of the total force F x as a function of wire separation d on 
wire 1 (a) and wire 2 (b) for a field incident from the right (— x direction) as in Fig. El Six different 
wavelengths are shown, and labeled in the legends. While the wires experience no net force in the 
y-direction, the forces in the x-direction vary greatly with wavelength and separation d. It is seen 
that depending on the wavelength, there are resonant peaks in the attractive forces between the 
two wires in the A range 410 — 470 nm. For the shorter wavelengths, A = 375 and 385 nm, there is 
relatively little response of the wire pair when compared to the longer wavelengths. The inset plots 
show the net induced y-component of the dipole moment, p y , corresponding to each force curve. 
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FIG. 6: The normalized ^-component of the force on wire 1 versus separation d. The force on the 
adjacent wire (wire 2) is identical except for a sign change. The incident field is directed downward 
as seen in the Fig. ^ The inset is the corresponding (normalized) dipole moment and the curves 
correlate with the legend and the same set of parameters are used as in the main plot. 

wavelengths, A = 375 and 385 nm where the F x curves remain negative and comparatively 
featureless over the entire d range. Clearly no resonant condition occurs for this shorter 
wavelength as is also evidenced in Fig. El These results are consistent with the two lower 
plots in Fig. 01 that show little gap field enhancement at A = 385 nm and a large gap field 
for A = 470 nm. The inset plots in Fig. El are similar in that the maxima in \F X \ show up as 
discernible peaks in \p y \, reflecting the charge redistribution. In general there is a redshift in 
the peak values and a sharp decline of \F X \ as d increases, indicating a corresponding shift 
in the resonance parameters of the mutually interacting wires. 

The normalized force, F x , on the left wire (wire 1) as a function of d for an incident 
field propagating downwards (in the — y direction) is shown in Figure |H1 This figure is 
analogous to Fig. As in Fig. [H the force on each wire is equal in magnitude and oppositely 
directed, yielding an attractive nature over the relatively broad wavelength range studied. 
The main resonant peaks are present at small separations and then F x monotonically declines 
towards zero with increased d. The wavelengths that give the largest F x occur in the 
range 385 — 410 nm. Outside of that range, increasing A has the effect of reducing the 
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force and narrowing its peaks. For sufficient separation distances, the force becomes weakly 
dependent on the wavelength. These characteristics reveal that there are optimal separations 
and wavelengths to utilize the plasmon driven force enhancement. The inset shows the 
normalized x-component of the dipole moment which is identical in each nanowire. The 
curves each follow a similar trend, with large separations maintaining a negative dipole 
moment, while for A = 375 nm the possibility exists to induce a positive net dipole moment 
for small separations. Thus for a particular value of d it is possible to have a vanishing 
dipole moment in the wire. This follows from regions within the wire structure that have 
oppositely oriented polarization vectors that sum to zero. 

In conclusion, we have studied the electromagnetic response of two long silver nanowires 
illuminated by a plane wave with the electric field perpendicular to the axis of the wires. 
The volume integral approach used here provides accurate solutions to Maxwell's equations 
over length scales that are much smaller than the wavelength of the driving field. We found 
that by varying the wavelength, separation distance, and angle of incidence of the electric 
field, the collective oscillation of the electrons induced appreciable forces and coupling in 
the nanowires. The results shown here can play a role in the fabrication and design aspects 
of optical spectroscopy and waveguiding over nanometer length scales. 
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